Diffusive Spreading of Chainlike Molecules on Surfaces 
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We study the diffusion and submonolayer spreading of 
chainlike molecules on surfaces. Using the fluctuating bond 
model we extract the collective and tracer diffusion coeffi- 
cients D c and Dt with a variety of methods. We show that 
D c (8) has unusual behavior as a function of the coverage 9. 
It first increases but after a maximum goes to zero as 8 — > 1. 
We show that the increase is due to entropic repulsion that 
leads to steep density profiles for spreading droplets seen in 
experiments. We also develop an analytic model for D c (9) 
which agrees well with the simulations. 
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The spreading dynamics of molecularly thin oil films on 
solid substrates has gained substantial interest recently 
This is not only due to its obvious significance 
for technical applications, but also to the richness of the 
structures (layering, fingering) in the observed film pro- 
files that sensitively depend on the molecular structure 
and interactions. As a typical feature, the development 
of molecularly thin so-called precursor films is observed 
0. The precursor film profile may be smooth, or ex- 
hibit pronounced molecular layering effects ("dynamical 
layering"). 

At the foremost tip of the film, the flow dynamics de- 
velops into surface diffusion of single molecules, i.e., a 2D 
gas. In some experiments it has been found that 

the transition from the dense submonolayer to the dilute 
film occurs much more rapidly than expected from simple 
molecular diffusion. The measured profiles are not Gaus- 
sian but assume a steeper shape that can be well fitted 
by a spherical cap || (see also Ref. 0]). An explanation 
for the late-time profiles has been suggested in terms 
of a transition of the fluid from a non-volatile 3D to a 
volatile 2D phase where the attractive intermolecular in- 
teractions are weaker |^,f|]. However, it is the aim of the 
present study to provide an alternative microscopic ex- 
planation for the steep density profiles observed in cases 
where desorption of the molecules can be neglected. By 
studying the diffusive properties of chainlikc molecules, 
we demonstrate that such profiles can be generated by 
entropic repulsion alone. 

To this end, we have performed a systematic study of 
diffusion and submonolayer spreading of athermal, flexi- 
ble chains. We concentrate on the coverage dependence 
of collective and tracer diffusion coefficients D c {6) and 



Dt(9), respectively. In addition to being important for 
submonolayer spreading dynamics, the diffusion of such 
molecules is of fundamental theoretical interest. Despite 
considerable experimental || and theoretical JtJ work on 
adatom diffusion on surfaces there are only few studies 
of more complicated molecules . 

We use the fluctuating bond (FB) model with Monte 
Carlo (MC) simulations to extract both D c (9) and D t {9) 
from the relevant correlation functions in equilibrium. 
These results are complemented by analysis of simu- 
lated density profiles together with the nonlinear diffu- 
sion equation. We show that while D t (9) is a decreasing 
function of 9 as expected, D c {9) displays more complex 
behavior. It first increases with 9, but after a maximum 
value goes to zero as 9 — ► 1. Wc show that it is this pos- 
itive slope of D c {9) that leads to steep precursor profiles 
in accordance with experiments gj. We develop a ther- 
modynamic theory for D c (9) which demonstrates that 
the initial increase is due to strong entropic repulsion, 
which is eventually overcome by the decreasing mobility 
of individual chains. The theory gives good agreement 
with our numerical results. 

The FB model [§ is a 2D lattice model of polymer 
chains, where each segment excludes four nearest and 
next-nearest neighbor sites on a square lattice. The bond 
lengths bg are allowed to vary between 2 < bg < Vl3 (in 
units of the lattice constant), where the upper limit pre- 
vents bonds from crossing each other. The stiffness of the 
chains is controlled through an angle dependent potential 

U$ = — Y17=i S^2i B-1 cos ( ( / ) ); where n is the number of 
chains, Npg is the number of segments in each chain and 
<fi the angle between two adjacent bonds. Dynamics is in- 
troduced in the model by Metropolis moves of single seg- 
ments, with a probability of acceptance min[e -A ' 7 */ T , 1], 
where AU^ is the energy difference. Only such moves are 
allowed that obey both the site exclusion and the bond 
length restrictions. One MC time step is defined as one 
attempt to move each monomer of every chain. We note 
that since no global translational modes of the chains are 
included, the rigid rod limit T — > is not well defined in 
the model go|. 

Simulations of diffusion were first done for T = 00 
using a 180 x 180 square lattice. We calculated Dt from 
the definition 
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where r*j(i) is the position of i th particle at time t. The 
results for Nfb = 6 are shown in the inset of Fig. 1. As 
expected, D t (8) is a strongly decreasing function of 8 due 
to the mutual blocking of the chains. 

The 

collective diffusion coefficient, however, shows strikingly 
different behavior. It can be defined through the Green- 
Kubo relation § D c = (l/2({Sn) 2 )) /* dt(J(0) ■ J(t)) 

where J(t) — Y17=i^i(t) 1S the total diffusion current, 
and {{5n) 2 ) the mean square fluctuation (in a finite area 
A). Since the definition of D c involves cross-particle 
velocity-velocity correlations, it samples the average col- 
lective density fluctuations of the chains instead of just 
single-chain properties. In Fig. 1 we show results for 
Nfb — 6 conveniently obtained using the temporal de- 
cay of the Fourier transformed density autocorrelation 
function S(k,t) = S^^c^ fn| where S is a con- 
stant. Great care was taken to ensure that the hydrody- 
namic limit J(| was reached. Initially, D c {8) increases up 
to 9 rs 0.7, after which it rapidly approaches zero. 

The initial rise of D c (9) is in agreement with conclu- 
sions drawn from the experimental studies [Q and MC 
simulations of density profiles Jl2],[l3| ■ We did additional 
studies simulating a circular 2D droplet spreading us- 
ing the FB model, and numerically solving the nonlinear 
diffusion equation d8(x,t)/dt — J^[D c (8)d9(x,t)/dx] by 
using a monotonically increasing tanh fitting function to 
D c (8) for 8 < 0.7. By matching the simulated profiles to 
the solutions of the equation, the form of D c (8) can be 
determined. We found it to be in perfect agreement with 
the MC results of Fig. 1. In Fig. 2 we show a comparison 
between experimental simulated, and numerically cal- 
culated density profiles in the submonolaycr regime. The 
agreement is excellent, and these profiles can be well fit- 
ted by a spherical cap shape with a Gaussian foot at 
lowest coverages j|. 

Next, we studied the effects of chain length and stiff- 
ness to collective diffusion. To study the full coverage 
range, we performed careful Boltzmann-Matano analy- 
sis of spreading profiles in a rectangular geometry. 
The results for Nfb = 6, T = oo were checked against 
the MC data, and very good agreement was found. In 
Fig. 1 we show these data for Nfb = 24 and 48, with 
T = oo. As expected, diffusion slows down but the qual- 
itative behavior remains the same. In Fig. 3 we show 
the effect of stiffness for Nfb = 6. Again, diffusion slows 
down but now also the maximum of D c {8) becomes less 
pronounced. 

To explain the unusual behavior of D c {8), we consider 
a simple thermodynamic theory for collective diffusion 
jjjj. The strong temperature dependence of D c (8) evi- 
dent in Fig. 2 suggests that besides the mobility of indi- 
vidual chains, entropy plays an important role here. Our 
starting point is the Nernst-Einstein equation §| , which 
relates collective diffusion to mobility m via 

D c (d) = me(^y (2) 



where fj, is the chemical potential. For chainlikc 
molecules, we approximate it by /i(8) = fiq + fcsT[ln(p) — 
ln(u;)] where p — 8/N is the number of molecules with 
N segments per unit area, and w is the number of con- 
figurations of the chain. This gives 

D c {8)=mk B T[l-8^^-]. (3) 

Next, both m and w must be estimated as a function of 
8. These quantities are model dependent. In general, the 
mobility is expected to decrease due to mutual blocking 
of the chains, and may be approximated by m w too[1 — 
c n {8)] ~ mo(l — 8 X ), where tuq is a constant and A is a 
scaling exponent for the probability of nearest neighbor 
occupation c„ . Using Flory theory in 2D it can be argued 
that A « 2.5 [^2|. For the entropy, in the case of a lattice 
polymer with N segments (N 3> 1) we can similarly write 
w « [1 + q(l — c n )] N , where q > 1 is a parameter. With 
these approximations, D c becomes 

where Do is the single chain diffusion constant which 
for single segment dynamics scales oc 1/N. Thus, we 
can write Dq = D s /N. In Fig. 4 we show the behavior 
of this expression as a function of coverage for various 
values of N. The curves are strikingly similar to the 
results obtained for the FB model, and the maximum in 
D c (8) becomes relatively more pronounced for increasing 
N Q. 

The theory above reveals that the underlying reason 
for the unusual behavior of D c is entropic repulsion, 
which at higher coverages is overcome by decreasing mo- 
bility of single chains. For the FB model with Nfb = 6, 
we have numerically extracted an effective pair interac- 
tion potential V e (r) []l6f which indeed shows a strong re- 
pulsion extending up to several lattice sites. We also 
find that all pair correlation functions G(r) for coverages 
8 < 0.7 and Nfb = 6,24 and 48 collapse to a single 
function which is given by g(x) = G(rd a /Np B ), where 
a w 0.38 and (3 w 0.55. Finally, it is interesting to note 
that the maximum of D c (8) occurs at coverages close to 
typical critical fractional coverages for continuum perco- 
lation of 2D objects Q. 

We can use the theory derived above to make a quan- 
titative comparison with the FB model. To this end, 
we have calculated the mobility and find that m 
D t (0)(l — 8 1A ) is a rather good approximation Jia]. To es- 
timate the entropy term for segments, we use the mean- 
field approximation w — w 2 W2~ 2 where W\ is the en- 
tropy arising from a segment at the end of the chain, and 
ui2 from one in the middle of the chain. For Nfb = 6, 
T = oo we have numerically determined w± (8) and Wi (8) , 
and find that they can be well fitted by a simple tanh 
function [^9|. In Fig. 1 the solid line indicates the re- 
sult for D c (8) as obtained from Eq. (3) where the only 
fitting parameter is N. We obtain very good agreement 
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with TV = 16 which is reasonable since the true entropy 
is underestimated by our approximation. 

To summarize, we have presented a combination ol nu- 
merical and analytic arguments to explain, how steep 
submonolayer density profiles observed in experiments 
H can be obtained solely from the entropic diffusion and 
spreading of chainlike molecules on surfaces. For ather- 
mal molecules, a strong entropic repulsion dominates at 
small and intermediate coverages, leading to the increase 
of D c (9). At high coverages, the chain mobility takes 
over and D c (9) approaches zero. The maximum of D c (6) 
is more pronounced for longer chains, but its position is 
relatively weakly dependent on the details of the model. 
This work also suggests that D c (8) can be used to obtain 
information about effective chain-chain interactions. 
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[15] It should be noted that Eq. (4) does not give the correct 
result for point-like particles, where D c = const. This 
limit can be built in by writing n(6) = u,o + kBT\n[6/ (1 — 
ff)\ which is exact for a Langmuir gas M, and m ~ D§{\ — 
6). 

[16] L. Verlet, Phys. Rev. 165, 201 (1968). 
[17] G. E. Pike and C. H. Seager, Phys. Rev. B 10, 1421 
(1974). 

[18] For the present system, this shows that setting D t ~ m 
H is a poor approximation and correlation effects are 
important. 

[19] For the FB model, fitting to 1 + q( 1 - 6 X ) gives unrealistic 
values for q and A, and the functional form is not well 
satisfied. 

FIG. 1. Collective diffusion coefficient vs. coverage for the 
FB model at T = oo. Points are MC data for Nfb = 6, 
and solid line is the theoretical fit with N = 16. Dashed 
and dash-dotted lines are Botzmann-Matano results for 
Nfb = 24 and 48. Inset shows the tracer diffusion coeffi- 
cient for Nfb = 6. The curves have been normalized by Di 
which is the diffusion coefficient for a single segment. 
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FIG. 2. Comparison between experimental density profiles 
of polydimethylsiloxine spreading on silver (circles), simula- 
tions of 2D circular droplets from the FB model with Nfb = 6 
(squares), and numerical solutions of the ID nonlinear diffu- 
sion equation (solid lines). 

FIG. 3. Influence of chain stiffness for collective diffusion 
for Nfb = 6. The data have been obtained from Boltz- 
mann-Matano analysis. 

FIG. 4. D c (8)/D s from Eq. (4), with A = 2.5 and q = 6. 
The solid lines from top to bottom (at = 0) are for 
N = 2, 6, 12, and 48. The maximum becomes more pro- 
nounced for large N, where Eq. (4) approaches a limiting 
form shown by the dashed line. 
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